Constraining dynamical dark energy with a divergence-free parametrization in the presence of 

spatial curvature and massive neutrinos 
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In this paper, we report the results of constraining the dynamical dark energy with a divergence-free param- 
eterization, w(z) = w + w a ( ln f^_ J - In 2), in the presence of spatial curvature and massive neutrinos, with the 
7-yr WMAP temperature and polarization data, the power spectrum of LRGs derived from SDSS DR7, the Type 
la supernova data from Union2 sample, and the new measurements of Hq from HST, by using a MCMC global 
fit method. Our focus is on the determinations of the spatial curvature, Q.^, and the total mass of neutrinos, 2 m v , 
in such a dynamical dark energy scenario, and the influence of these factors to the constraints on the dark energy 
parameters, wo and w a . We show that Q.^ and J] ni v can be well constrained in this model; the 95% CL limits 
are: -0.0153 < Sl k < 0.0167 and Y, m v < 0.56 eV. Comparing to the case in a flat universe, we find that the 
error in w is amplified by 25.51%, and the error in w a is amplified by 0.14%; comparing to the case with a zero 
neutrino mass, we find that the error in wo is amplified by 12.24%, and the error in w a is amplified by 1.63%. 

PACS numbers: 95.36.+X, 98.80.Es, 98.8().-k 
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I. INTRODUCTION 

Dark energy is one of the most important themes in physics 
today. However, we do not know much about dark energy 
due to the accuracy of current data. Though the current ob- 
servational data are consistent with a cosmological constant, 
the possibility that dark energy is dynamical is still not ex- 
cluded by the data and has been attracting wide attentions in 
the cosmology and theoretical physics communities. 

In order to detect the dynamics of dark energy, one usually 
has to parameterize the equation-of-state parameter (EOS), w, 
empirically, using two or more free parameters. Among all 
the parametrization forms of EOS, the Chevallier-Polarski- 
Linder (CPL) parametrization [1], w(a) — wq + w a (l - a), 
where wq and w„ are parameters and a is the scale factor of 
the universe, is the most widely used one and has been ex- 
plored extensively. However, as pointed out in Ref. [2], the 
CPL description will lead to unrealistic behavior in the future 
evolution, i.e., \w(z)\ grows rapidly and eventually encounters 
divergence as the redshift z approaches - 1 . In order to keep 
the advantage of the CPL parametrization and avoid its draw- 
back at the same time, it is believed that a parametrization that 
is free of divergence both in the past and future evolutions is 
necessary. 

In Ref. [2], Ma and Zhang proposed the following hybrid 
form of logarithm and CPL parametrizations: 



/ln(2+z) 
w(z) = wo + w a I - - In 2 



(1) 



This novel parametrization has well behaved, bounded evolu- 
tion for both high redshifts and negative redshifts. In particu- 
lar, for the limiting situation, z — > -1, a finite value for EOS 
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can be obtained, w — wq + w a (l - In 2). At low redshifts, 
this parametrization form reduces to the linear one, w(z) ~ 
wo + w a z, where w a = -(ln2)w fl . Of course, one can also re- 
cast it at low redshifts as the CPL form, w(z) ~ wo+w a z/(l+z), 
where w a = (1/2 - ln2)w n . Therefore, it is clear to see that 
this parametrization exhibits well-behaved feature for the dy- 
namical evolution of dark energy. Without question, such a 
two-parameter form of EOS can genuinely cover scalar-field 
models as well as other theoretical scenarios. In Ref. [3], this 
parametrization was used to explore the ultimate fate of the 
universe. 

The parametrization (1) has been explored deeply by us- 
ing the current data. In particular, in Ref. [4], we have an- 
alyzed the detection of dynamics of dark energy with the 
parametrization (1) by performing a Markov Chain Monte 
Carlo (MCMC) global fitting method. Since dark energy pa- 
rameters are tightly correlated with the spatial curvature 
and neutrino mass £ m v , in this paper we will deeply analyze 
the influences of these factors on the detection of dynamics of 
dark energy. 

The paper is organized as follows: In Sec. II we will in- 
troduce our global fitting procedure and the data we used for 
analysis; the results are presented in Sec. Ill, and our conclu- 
sion is given in Sec. IV. 



II. GLOBAL FITTING PROCEDURE AND DATA 

We have modified MCMC package CosmoMC [5] to per- 
form a global fitting analysis for the dynamical dark energy 
parameterized above. As we know, within the dynamical dark 
energy models, the perturbations of dark energy are impor- 
tant [6-8] for the data fitting analysis. For quintessencelike or 
phantomlike models, whose w does not cross the cosmologi- 
cal constant boundary, the perturbation of dark energy is well 
defined. However, when w crosses -1, which is described 
by a quintom dark energy model [9], one is encountered with 
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the divergence problem for perturbations of dark energy at 
w = - 1 . For avoiding such kind of divergence problem, in this 
paper we use the method provided in Refs. [8, 10] to treat the 
dark energy perturbations consistently in the whole parameter 
space in the numerical calculations. 

Our most general parameter space vector is: 

P = (a)b, co c , 0, x, wo, w a , Q.k, ^ m v , n s ,A s , c 2 s ), (2) 

where a>h = £lbh 2 and to c = Q. c h 2 , with and Q. c the physi- 
cal baryon and cold dark matter densities relative to the critical 
density, Q. k is the spatial curvature satisfying + Cl,„ + fide 
I . Yt m v is the total mass of neutrinos, is the ratio (multi- 
plied by 100) of the sound horizon to the angular diameter 
distance at decoupling, x is the optical depth to re-ionization, 
wo and w a are the parameters of dark energy EOS given by 
Eq. (1), A s and n s are the amplitude and the spectral index 
of the primordial scalar perturbation power spectrum, and c s 
is the sound speed of dark energy. For the pivot scale we set 
k s o = 0.05Mpc~'. Note that we have assumed purely adiabatic 
initial conditions. 

Note that the sound speed of dark energy is fixed in our 
analysis. In the framework of the linear perturbation theory, 
besides the EOS of dark energy, the dark energy perturbations 
can also be characterized by the sound speed, c 2 s = 6pd e /Spd e - 
The sound speed of dark energy might affect the evolution of 
perturbations, and might leave signatures on the CMB power 
spectrum [11]. However, it has been shown that the con- 
straints on the dark energy sound speed c 2 in dynamical dark 
energy models are still very weak, since the current observa- 
tional data are still not accurate enough [12]. Therefore, in 
our analysis, we have treated the dark energy as a scalar-field 
model (multi-fields or single field with high derivative) and 
set c 2 to be 1. Of course, one can also take c 2 as a parameter, 
but the fit results would not be affected by this treatment [12]. 

In the computation of the cosmic microwave background 
anisotropy (CMB), we include the 7-year WMAP tempera- 
ture and polarization power spectra [13] with the routine for 
computing the likelihood supplied by the WMAP team [14]. 
For the large scale structure (LSS) information, we use the 
power spectrum of luminous red galaxies (LRGs) measured 
from the SDSS DR7 [15]. The supernova (SN) data we use 
are the recently released "Union2" sample of 557 data [16]; 
note that the systematic errors are included in our analysis. In 
the calculation of the likelihood from SN we marginalize over 
the relevant nuisance parameter [17]. 

Furthermore, we make use of the Hubble Space Tele- 
scope (HST) measurement of the Hubble constant Ho = 
lOO/i km s _1 Mpc~' by a Gaussian likelihood function cen- 
tered around h = 0.738 and with a standard deviation cr 
0.024 [16]. 



III. NUMERICAL RESULTS 

In this section we shall present the results of the global fit- 
ting analysis. In particular, we shall test the influences on 
probing the dynamics of dark energy when spatial curvature 



and massive neutrinos are considered in the data analysis. 
Since we know that both the spatial curvature parameter Q.^ 
and the total mass of neutrinos 2 m v are degenerate with the 
equation of state of dark energy w, we are very much inter- 
ested in seeing how these two factors affect the constraints on 
the parameters wo and w a . Moreover, we also want to see how 
the current data constrain and £ m v within the framework 
of dynamical dark energy with parametrization (1). 



A. Spatial curvature and EOS 

Dark energy parameters and D.^ are correlated via the cos- 
mological distance information. From the observation of 
CMB, the curvature of the observable universe can be de- 
termined by the position of first acoustic peak of CMB tem- 
perature power spectrum precisely. However, Q. k is tightly 
degenerated with f2,„, and such degeneracy can be, in certain 
level, broken by taking into account the data of large scale 
structure and supernovae. 

In constraining Q.^, one usually uses a combination of dis- 
tance measurements from baryon acoustic oscillations (BAO) 
and CMB. In this paper, however, we use the power spectrum 
of LRGs measured from SDSS DR7, instead of BAO, to con- 
strain Qfc. Note that BAO and the LRG power spectrum cannot 
be treated as independent data sets because a part of the mea- 
surement of BAO comes from LRGs as well. Since the LRG 
power spectrum is a powerful probe of the total mass of neu- 
trinos, 2 m v, in this paper we uniformly use the LRG power 
spectrum, instead of BAO, to constrain and £ m v . In the 
following we shall consider the case of massless neutrinos and 
see the influence of fi^. 

In Fig. 1, we present the two dimensional cross correlation 
and one dimensional probability distribution of f2<., wo, w n , 
D. m and H by fitting with 7-yr WMAP, LRG, SN and H data. 
From the panels of Q*-wo and Q.k-w a , we find that the corre- 
lation between and wo is positive, but we do not observe 
much correlation between and w fl . 

The 68% confidence level (CL) constraint on is 0.0003 ± 
0.0079, from which we can see that our universe is very close 
to flatness. The 95% CL limit is: -0.0153 < Q A < 0.0167. So, 
we find that in the dynamical dark energy model (1) the spatial 
curvature is well constrained by WMAP+LRG+SN+// . It is 
of interest to make a comparison to the cases of ACDM (w = 
-1) model and constant w model. Komatsu et al. gave the 
limit on £l k for the case of w = -1, Q. k = -0.0023+° ™^ (68% 
CL), from WMAP+BAO+Ho, and for the case of constant 
w, Qjt = -0.0057^^6 (68% CL), from WMAP+BAO+SN, 
where SN is the Constitution sample. These results are in 
good agreement with our limit. 

How does the spatial curvature parameter, Q. k , affect the 
constraint results of (wo, w a )l We find that, in a non-flat uni- 
verse, the results are: 

wo = -0.922 + 0.123 and w a = 1.651 + 1.470 (68% CL). 

Thus, even when w is allowed to depend on time, the cur- 
rent data are still consistent with a cosmological constant. 
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FIG. 1: Constraints on the cosmological parameters Qj, Wo, w , fl, 
and Hg m tne dynamical dark energy scenario (1). 



However, a large range of values of (wo, w a ) are still al- 
lowed by the data. Comparing to the case in a flat universe: 
w = -0.921 + 0.098 and w a = 1.905 + 1.468 (68% CL), 
we find that the best-fit values of wo and w a are shifted by 
0.1% and 13.3%, respectively; the error in wo is amplified by 
25.5 1 %, and the error in w a is almost the same, only amplified 
by 0.14%. 

One may be curious to know why the error on wo is signif- 
icantly amplified, but that for w a is scarcely affected by the 
presence of O^. Actually, this is fairly easy to be understood. 
From the fitting results, it is clear to see that wo is constrained 
more tightly than w fl , i.e., wo is more sensitive to the data, 
comparing to w a . The current constraint on w a is very weak, 
only at the level of ~ 80%, while for wo the constraint is much 
stronger, at the level of ~ 10%. Thus, once Q* is involved in 
the analysis, its influence on wo is, obviously, much greater 
than that on w a . 



B. Neutrino mass and EOS 

Weighting the neutrino mass is one of the most important 
challenges in modern physics. Currently the neutrino oscilla- 
tion experiments, such as atmospheric neutrinos experiments 
[18] and solar neutrinos experiments [19], have confirmed 
that the neutrinos are massive, but give no hint on their ab- 
solute mass scale. Cosmology can provide crucial comple- 
mentary information on absolute neutrino masses, because 
massive neutrinos leave imprints on the cosmological obser- 
vations, such as the Hubble diagram, CMB temperature power 
spectrum and LSS matter power spectrum [20]. 




FIG. 2: Constraints on the cosmological parameters vfo, w a , fl,„, Ho 
and Yi m v in the dynamical dark energy scenario (1). 



In Ref. [13], Komatsu et al. have constrained the to- 
tal mass of neutrinos, 2 m v = 94eV(Q v /i 2 ), from the 7-yr 
WMAP data combined with the distance information. For 
a flat ACDM model, i.e., w = -1 and Q.^ = 0, they found 
that the WMAP+BAO+# limit is £/n„ < 0.58 eV (95% 
CL). The limit improves when information on the growth of 
structure is added. For example, they found that, when the 
BAO is replaced by the power spectrum of LRGs, the combi- 
nation WMAP+LRG+i/o gives £ m v < 0.44 eV (95% CL) 
for w = -1. Thus, we can see that the power spectrum 
of LRGs plays an important role in constraining the neu- 
trino mass. In the following we shall use the combination of 
WMAP+LRG+SN+// to constrain the neutrino mass, £ m v , 
as well as the parameters of dynamical dark energy, wo and 
w a , in a flat universe. 

In Fig. 2, we present the two dimensional cross cor- 
relation and one dimensional probability distribution of 
2 m v , wo, w fl , Q,„ and Hq, from the data combination 
WMAP+LRG+SN+// - From the panels of ^My-Wo and 
2 my-w a , we find that there is an anti-correlation between wq 
and 2 "iy, but no significant correlation between w a and £ m v 
is observed. 

The global fitting gives the constraint on the neutrino mass 
in our dynamical dark energy scenario: 2 m v < 0.56 eV (95% 
CL). So, it is interesting to find that in our dynamical dark 
energy scenario (1) the total mass of neutrinos can be well 
constrained by the data combination WMAP+LRG+SN+// . 
To make a comparison, we show the results of the constant w 
model given by Komatsu et al. [13]: 2 m v < 0.71 eV (95% 
CL) from WMAP+LRG+// , and 2 m v < 0.91 eV (95% CL) 
from WMAP+BAO+SN (where SN is the Constitution data). 
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In our case, even though there is one more parameter, w a , the 
constraint on £ m v is much tighter. This shows that the data 
combination we used, WMAP+LRG+SN+// , is fairly good 
at constraining the neutrino mass. To show this more clearly, 
we also use our data combination to make a global fitting cal- 
culation for the case of the constant w model, and the con- 
straint result we obtained is £ tn v < 0.45 eV (95% CL). 

How does the neutrino mass, £ m v, affect the constraint re- 
sults of (wo, w a )l We find that, with the non-zero neutrino 
mass, the results are: 

w = -0.972 + 0.110 and w a = 2.038 + 1.492 (68% CL). 

So, still, though the current data are consistent with a cosmo- 
logical constant, a large range of values of (wo, w fl ) are still 
allowed by the data. Comparing to the case with a zero neu- 
trino mass, wo = -0.921 +0.098 and w a = 1 .905 + 1 .468 (68% 
CL), we find that the best-fit value of wo is shifted by 5.5%, 
and the best-fit value of w„ is shifted by 7.0%; the error in wo 
is amplified by 12.24%, and the error in w a is amplified by 
1.63%. Obviously, the reason that the influence on the error 
of wo by the presence of £ m v is much severer than that of w„ 
is as the same as the case of the spatial curvature: wo is much 
more sensitive to the data than w fl . 



IV. CONCLUSION 

With the 7-yr WMAP temperature and polarization data, 
the power spectrum of LRGs derived from SDSS DR7, the 
Type la supernova data from Union2 sample, and the new 
measurements of Ho from HST, we have tested the dynamical 
dark energy parametrization, w(z) = wo + w a ( in f^ - In 2), 
in the presence of spatial curvature and massive neutrinos, by 
using a MCMC global fit method. For a time-dependent equa- 
tion of state, one must be careful about the treatment of per- 
turbations in dark energy when w crosses -1. We used the 
method provided in Refs. [8, 10] to treat the dark energy per- 
turbations consistently in the whole parameter space in our 
numerical calculations. The sound speed of dark energy, c 2 s , 
is fixed in our calculation to be 1; note that the value of c 2 s 
is actually insensitive to our final fit results. Our focus was 
put on the determinations of the spatial curvature of the ob- 
servable universe, Q^, and the total mass of neutrinos, £ m v , 
in our dynamical dark energy scenario, and the influence of 



these factors to the constraints on the dark energy parameters, 
wo and w a . 

The 95% CL limit on the spatial curvature is -0.0153 < 
Qk < 0.0167. Thus, we find that in our dynamical dark 
energy model the spatial curvature is well constrained by 
WMAP+LRG+SN+7/o- In a curved universe, the constraint 
results of wo and w a are wo = -0.922 + 0.123 and w a = 
1.651 + 1.470 (68% CL). Comparing to the case in a flat uni- 
verse, w = -0.921 + 0.098 and w a = 1.905 + 1.468 (68% 
CL), we find that the error in wo is amplified by 25.51%, and 
the error in w a is almost the same, only amplified by 0.14%. 
Our global fit gives the constraint on the neutrino mass in 
our dynamical dark energy scenario, 2'«v < 0.56 eV (95% 
CL). So, we find that in our dynamical dark energy scenario 
the total mass of neutrinos can be well constrained by the 
data combination WMAP+LRG+SN+H . With the non-zero 
neutrino mass, the constraint results of wo and w fl are wo = 
-0.972 ± 0.110 and w a = 2.038 + 1.492 (68% CL). Compar- 
ing to the case with a zero neutrino mass, wo = -0.921 ±0.098 
and w a = 1.905 ± 1.468 (68% CL), we find that the error in 
wo is amplified by 12.24%, and the error in w a is amplified by 
1.63%. 

We believe that it is fairly important to use some 
divergence-free parametrization to probe the dynamical evo- 
lution of dark energy. We have shown that the form (1) is a 
good proposal, and it has been proven to be very successful in 
exploring the properties of dark energy. We suggest that this 
parametrization should be further investigated. 
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